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Abstract 


A two-dimensional model comprising fuel channel, anode, cathode and electrolyte layers for anode-supported micro-tubular solid oxide fuel cell 
(SOFC), in which momentum, mass and charge transport are considered, has been developed. By using the model, tubular cells operating under 
three different modes of current collection, including inlet current collector (IC), outlet current collector (OC) and both inlet and outlet collector 
(BC), are proposed and simulated. The transport phenomena inside the cell, including gas flow behavior, species concentration, overpotential, 
current density and current path, are analyzed and discussed. The results depict that the model can well simulate the diagonal current path in 
the anode. The current collecting efficiency as a function of tube length is obtained. Among the three proposed modes, the BC mode is the most 


effective mode for a micro-tubular SOFC, and the IC mode generates the largest current density variation at z-direction. 


© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


Solid oxide fuel cells (SOFCs) are one of the most promising 
power generation technologies because of their high efficiency 
and clean generation of electric power [1,2]. Tubular SOFCs 
have several advantages compared with planar ones. For exam- 
ple, the sealing area is much smaller in a tubular cell than in a 
planar cell [3]. Fuel and air distribution manifolds can be more 
facile in tubular stacks. Thermal stresses can be less troubled 
in tubular systems because each membrane-electrode assembly 
is less rigidly restricted. Recent years, anode supported tubular 
SOFC has been rapidly developed [4—8]. The much longer cur- 
rent path in a tubular SOFC than in a planar cell becomes one 
of the most concerned issues. Especially for an anode supported 
micro-tubular SOFC, current at the anode side can only be drawn 
from anode ends. The current collecting mode is very important 
in tubular cells for achieving high cell efficiency. 


Abbreviations: BC, both inlet and outlet collector; EMF, electromotive force; 
IC, inlet current collector; OC, outlet current collector; SOFC, solid oxide fuel 
cell; TC, total area current collector 
* Corresponding author. Tel.: +86 411 84379049; fax: +86 411 84379049. 
E-mail address: mjcheng @ dicp.ac.cn (M. Cheng). 
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Numerical calculations are a very useful tool in the SOFC 
research because the detailed measurement on the phenomena 
in SOFCs is difficult [9-11]. Up to now, most cell performance 
simulations have been carried out using the circuits of an equiv- 
alent circuit of a cell portion plus a resistance portion, namely, 
a ladder circuit [10-12]. The current-flow directions obtained 
from these models are inconsistent with the facts in a real cell. 
When there is a large difference of the electromotive force (EMF) 
and a long current path in a cell, the current may flow diagonally 
in the electrolyte and electrodes, and the ladder equivalent cir- 
cuit model cannot simulate the diagonal current in the electrolyte 
and electrodes precisely [12]. In order to deal with the diago- 
nal current in the electrolyte and electrodes, it is necessary to 
develop a new model for simulating the current path in the cell 
components. 

In this paper, a two-dimensional model comprising fuel 
channel, anode, cathode and electrolyte layers for anode- 
supported micro-tubular solid oxide fuel cell will be built. The 
calculations of electro-chemical, fluid phenomena and current 
path inside the single tubular cell are correlated to each other. 
We will investigate the current density at the anode/electrolyte 
interface and the current collecting efficiency under three dif- 
ferent current collecting modes of inlet current collector (IC), 
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Nomenclature 

c concentration (mol m~?) 

D diffusion coefficient (m? s7!) 

Eo standard potential (V) 

F Faraday’s constant (96485 C mol`!) 
i current density (A m~?) 

ig exchange current density (A m~?) 
I current (A) 

K permeability (m7) 

Ne electron transferred per reacting 
p pressure (atm) 

Rg gas constant (J mol`! K7!) 

T temperature (K) 

u velocity (m s7!) 


Veell cell voltage (V) 


Greek letters 


B transfer coefficient 

E porosity 

(0) potential (V) 

n overpotential (V) 

u dynamic viscosity (m? s7!) 
p density (kg m~?) 

o conductivity (S m7!) 
T tortuosity factor 
Subscripts 

a anode 

act activation 

c cathode 

diff diffusion 

ele electrolyte 

K Knudsen diffusion 
m ordinary diffusion 
rev reversible 

0 standard 


outlet current collector (OC) and both inlet and outlet collector 
(BC). 


2. Model description 
2.1. Basic feature of the cell and model geometry 


For a micro-tubular cell, current collector at anode side is only 
attached to the end of the anode tube. Hz humidified at room tem- 
perature is introduced as a fuel at the flow rate of 30cm? min™!. 
Oxygen concentration at the cathode/air interface is assumed to 
be the same as in air. Using this model, the cell performances 
with different tube length are calculated. The effects of three dif- 
ferent current collecting modes are analyzed. The dimensions of 
the cell are listed in Table 1. 

The anode supported tubular cell consists of YSZ electrolyte, 
Ni-YSZ anode and LSM-YSZ cathode. Because the micro- 


tubular cell has a symmetrical structure, it is assumed that the 


Table 1 
Input parameters to model [12,22] 
Anode conductivity (S m7!) Oa = 95x10 x exp Z =~) 
Cathode conductivity (S m7!) Oc = ae x exp (= zay 
Pap = 4 ( — +) 
Electrolyte conductivity (S m~*) Oe = 3.34 x 10° x exp 
Tubular interior diameter (mm) 1.49 
Anode thickness (mm) 0.72 
Cathode thickness (um) 50 
Electrolyte thickness (wm) 50 
Transfer coefficient (8) 0.5 
Electron transferred per reacting 1 
(ne) 
Anode exchange current density 5300 
(ipa) (A m°) 
Cathode exchange current density 2000 
Cios) (Am™”) 
K (m2) le—12 
um (Pas) 21.63e—7 
umo (Pas) 40.45e—7 
Do, c (m? s7!) 7.588e—6 
Dy; c (m? s7!) 7.588e—6 
Dehannel,m m? s7!) 8.506e—4 
Dp; a (m? s7!) 3.14e—5 
Dp,0,a (m? s7!) 1.39e—5 


cell properties are uniform in each angular section. Thus, one 
geometrical dimension is neglected, and the 2D modeling is con- 
sidered for the entire cell performance, as shown in Fig. 1. The 
current collector at the anode side is attached to one or two tube 
ends, whereas the current collector at the cathode side was con- 
nected to the whole cathode. The anode tube length and current 
collecting modes are crucial for the reduction of efficiency loss 
[8]. 

Three modes of current collection, including inlet current 
collector (IC), outlet current collector (OC), both inlet and outlet 
collector (BC) are proposed and calculated. In order to compare 
current collecting efficiency under the three different current 
collecting modes, the cell performance with all anode interior 
surfaces for current collection (TC) was also calculated as a 
reference. The current-collecting efficiency is defined as the cell 


Electrolyte 
Current collector 


Cathode 


Fig. 1. Schematic diagram of a micro-tubular geometry in 2D. 
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performance ratio of the proposed current collecting mode to the 
TC mode. 


2.2. Model assumptions 


e The model is based on steady state, using humidified H2 as 
the fuel. 

e Gas flows in the channels are laminar. 

e The reactant gas mixtures are approximated as ideal gas and 
incompressible. The fuel cell operates with 100% current 
efficiency. 

e The model is assumed to be isothermal, and the cell run at 
800°C. 


2.3. Electrochemical model 


The reactions to deal with the SOFC model are those com- 
monly reported in the literature [13,14]. The oxygen anions from 
oxygen reduction at cathode pass through the electrolyte and 
react with hydrogen at the anode-—electrolyte interface: 


102 + 2e7 > O77 (1) 
H +0% > H20 + 2e7. (2) 


Activation polarization losses occur from slow electron 
release or capture steps at the electrode—electrolyte interfaces. 
The Butler—Volmer kinetic equations are generally adopted for 
expressing the local current density distribution [15]. 


: ; Ne Fact Nef Nact 
ia = lao 9 exp | £ ) exp | (l=) | \ 
p { ( RoT RoT 
(3a) 


A . NeF act NeF act 
ic = ic,0 {exp (s ) exp | d — p) | \ 
RT RT 


ia. is 5300 A m~? for a Ni-YSZ anode, ic. is 2000 A m~? for a 
LSM-YSZ cathode. The activation loss can be calculated accord- 
ing to the following equation: 


(3b) 


Nact = Prev — |Pelectrode — Pele| (4) 


where rey is the reversible potential between the electrode 
and the electrolyte. @rey can be obtained from Nernst equa- 
tion. Since the electrochemical active layer is usually very thin 
compared with the electrode, the current generation is consid- 
ered to take place only at the electrode—electrolyte interfaces 
[14]. 

In expression (3a), iao, ico and £ are the exchange cur- 
rent density and the transfer coefficient (A m2). The transfer 
coefficient is considered to be the fraction of the change in 
polarization that results in a change in the reaction rate con- 
stant. The exchange current density is the forward and reverse 
electrode reaction rate at the equilibrium potential, and a high 
exchange current density means that a high electrochemical 


reaction rate and good fuel cell performance can be expected 
[16]. 
Electronic charge at the anode: 


—V(aV ba) = 0 (Sa) 

anode/electrolyte interface boundary: 

—o,Vdba-n = —ig (5b) 

anode current collector boundary: 

anode = 9 (5c) 

other anode boundaries: 

—0,Voa:n = 0 (5d) 
Ionic charge at the electrolyte: 

—V (Cele V dele) = 0 (6a) 


anode/electrolyte interface boundary: 


— Oele V Pele ` N = ia (6b) 
cathode/electrolyte interface boundary: 
—0ele V Pele : n = —ic (6c) 


other electrolyte boundaries: 
—OcleV Qele : n = 0 (6d) 
Electronic charge at the cathode: 
—V(o. Vdc) = 0 (7a) 
cathode/electrolyte interface boundary: 
—o.Vebe-n = ic (7b) 
cathode current collector boundary: 
dc = Veel (7c) 
other cathode boundaries: 
—o.Voe-n = 0. (7d) 
2.4. Momentum and mass transport equations 
To calculate the velocity distribution in the gas channel, the 


steady-state incompressible Navier-Stokes equation, Eq. (8a) 
and Eq. (8b) are solved numerically [17]: 


puVu) — Vu(Vu + (Vu)T) + VP =0 (8a) 
Vu =0 (8b) 


The inlet fuel is composed of H20 and H? (0.97:0.03 in molar 
fractions), flowing at a volume velocity of 30 sccm at room tem- 
perature. The density (p) is calculated from the ideal gas law, 
and the mean molar dynamic viscosity (jz) is obtained from the 
product of individual viscosity and molar fraction. An initial gas 
volume velocity of 30 sccm and a pressure (Po) of 101325 Pa are 
applied to the gas inlet and outlet boundaries Eqs. (8c) and (8d). 
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At the interfaces of electrolyte and electrode, no slip conditions 
are employed Eq. (8e). 


u-n = uo (8c) 
n(—Vu(Vu + (Vu)T) + VP) = —nPo (8d) 
u=0 (8e) 


Darcy equation describes the flow in the porous structure well 
far away from the walls, but it cannot model a no-slip condition 
near a wall. Thus, Brinkman equation is used for the porous 
media [18,19]. 


VP = uV°u + ral (9) 


where K is the permeability of the medium considered. 

Mass diffusion coefficients are required whenever species 
transport equations in multi-component flows are solved. Dif- 
fusion in porous medium is usually described by molecular 
diffusion or Knudsen diffusion. Knudsen diffusion occurs when 
the pore diameter is smaller than the mean free path of the gas 
molecules. Molecular diffusion occurs when the pore diame- 
ter is larger than the mean free path of the gas molecules. For 
an SOFC, both Knudsen and ordinary molecular diffusion pro- 
cess have to be considered [20]. The overall effective diffusion 
coefficient is given by [15] 


A eee a (10) 
‘ t\ Dim Dk 


where £, t, Dim and D x denotes porosity of porous 
medium, tortuosity factor, the ordinary and Knudsen diffusion 
coefficients[21], respectively. Finally, the steady-state diffusion 
and convection equation of Eq. (11) is solved from the velocity 
distribution to obtain the concentration distribution of reactants 
and products in the gas channel and the electrodes. 


V(—D;Vcj + ciu) = 0 (11) 


anode/electrolyte interface boundary: 


(—Dyy,aV CH + CHU) n= - 34 (12a) 
(Daa Voro + cmon) n = 3% (12b) 
anode channel inlet: 

CH) = CH3,0 (12c) 
CH0 = CH0,0 (12d) 
anode channel outlet: 

(—Dehannel, mV Ci + ciu) +n = ciju -n (12e) 
cathode/electrolyte interface boundary: 

(—Do,,cVco, + cou): n = -5 (13a) 


O2 concentration at the cathode/air interface is equal to the O2 
concentration in air as assumed at Section 2.1. 


CO. = CO>,0 (13b) 
other boundaries: 
(—D;Vcj + cju)-n = 0. (13c) 


2.5. Numerical implementation 


The calculations are solved using the finite element com- 
mercial software COMSOL MULTIPHSICS®, Version 3.3a. 
Because the primary objective in this modeling is to compare the 
cell performance under the different current collecting modes 
and the current density distribution, the mathematical simula- 
tion was performed at 0.7 V. The outputs from the model are the 
distributions of velocity, current density and species concentra- 
tions. 


3. Results and discussion 
3.1. Flow behavior and distribution of species 


From the simulation, it is evident that velocity decreases 
rapidly in the porous anode. Fig. 2 shows the calculated pro- 
file of the velocity for the fuel flow in the channel and the anode. 
Because the constant temperature and the incompressible gases 
are assumed, the velocity profile changes slightly along the z- 
direction. The velocity close to fuel channel wall is obviously 
falling due to the viscosity stress. In this simulation, the mixed 
gases of H2/H20 with a channel center velocity of about 2 m s7! 


Max: 2.064 
0.02 u[m/s] ax 
2 
1.8 
0.016 
1.6 
1.4 
0.012 
_ 1.2 
E£ 
N 1 
0.008 
0.8 
0.6 
0.004 
0.4 
0.2 
0 
r 0 
Min:0 


Fig. 2. Surface plot of velocity profile for fuel in anode and fuel channel for a 
2cm long SOFC (current collector at inlet). 
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(b) 
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(c) 


02 z: H2 Max: 0.970 os 
0.02 Ma Ns 0.02 0.02 H2O Max: 0.342 
0.95 
0.209 0.3 
0.016 0.016 0.016 
0.208 0.9 
0.25 
0.012 0.207 0.012 0.85 0.012 
z g g 0.2 
Š 0.206 N N 
0.8 
0.008 0.008 0.008 
0.205 0.15 
0.75 
0.204 
0.004 0.004 0.004 0.1 
0.203 0.7 
0.05 
0 0.202 0 0 
r Min: 0.202 r Min: 0.658 r Min: 0.0300 


Fig. 3. Surface plots of O2 (a), Hz (b) and H20 (c) and molar fraction distribution for a 2 cm long SOFC (only current collector at inlet). 


are imported from the SOFC fuel channel inlet at 800°C. The 
velocity of fuel flow in anode is reduced to about 1076 of that 
in the channel. 

Fig. 3(a) shows the distribution of oxygen molar fraction in 
the porous cathode. Due to the oxygen reduction reaction, the 
oxygen molar fraction decreases from the cathode/air interface 
to cathode/electrolyte interface. The distributions of hydrogen 
and water vapor molar fractions in the anode channel and porous 
anode are shown in Fig. 3(b) and (c). Since the hydrogen con- 
sumption occurs at the anode/electrolyte interface, the hydrogen 
molar fraction decreases along the fuel-flow direction, whereas 
the water vapor molar fraction increases along the z-direction. 
The variations of hydrogen and water vapor molar fractions in 
the porous anode along r-direction are very sharp but very slow 
in fuel channel. These results show that species transport resis- 
tance of porous anode can lower the fuel gas diffusion velocity in 
porous anode and reduce the cell performance. In one word, the 
overall cell performance of anode-supported tubular SOFC can 
be limited by the fuel gas diffusion velocity in the thick porous 
anode. 

Fig. 4 shows the variation of H2 molar fraction along the 
anode/electrolyte interface. H2 molar fraction decreases along 
the fuel-flow direction under all the four current collecting 
modes. The H2 molar fraction at inlet is smaller under the IC 
mode than under the OC mode. The variation of H? molar frac- 
tion along the z-direction is the mildest under the IC mode. Thus, 
the electrochemical reaction at fuel inlet is faster under the IC 
mode than under the OC mode. It is also shown in Fig. 4 that the 
reaction at the outlet is relatively slow under the IC mode, and 
H2 molar fraction at outlet is the highest among four current- 
collecting modes. The variation of Hz molar fraction is much 
sharper under the OC mode than under the IC mode. The main 


reason is that the inlet electrochemical reaction is milder under 
the OC mode than under the other modes, and less H2 is con- 
sumed. However, the reaction at outlet is relatively faster under 
the OC mode than under the other modes, and more Ho is con- 
sumed at outlet. The Hz molar fractions under both the BC and 
TC modes are the same as that under the IC mode at inlet, and 
decreases from inlet to outlet with a rate similar to that under the 
OC mode, indicating a higher H? depleting rate under the two 
modes. It can be seen in Fig. 4 that the variation of H20 molar 
fraction along the anode/electrolyte interface is opposite to that 
of H2. 

Fig. 5 shows the changes of O molar fraction along the 
cathode/electrolyte interface. Since O2 molar fraction at the 
cathode/air interface is assumed to be equal to that in the air 
(0.21), the amount of O2 consumed in electrode can obviously 
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Fig. 4. Hz and H20 molar fractions along the anode/electrolyte interface for a 
2cm long SOFC. 
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Fig. 5. O2 molar fraction along the cathode/electrolyte interface for a 2 cm long 
SOFC. 


be seen. Under the IC, BC and TC modes, the consumed Oz frac- 
tions at the inlet are very close and higher than that under the 
OC mode, and the O2 molar fraction at electrolyte/cathode inter- 
faces increases from inlet to outlet. However, O2 molar fraction 
exhibits a parabola-like curve under the OC mode, indicating 
that the O2 consumption is quite even along the cell. 


3.2. Electrochemical features and distribution of current 
field 


Fig. 6 shows the anode, cathode electric potential, electrolyte 
potential distributions and the current direction in the electrolyte. 
Because the potential of the anode current collector is assumed 
to be zero, the electric potential in the porous anode is nega- 


(a) (b) 
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0.7 0.88 
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Min: 0.700 


tive and gradually decreases from the current collector to the 
anode/electrolyte interface. Due to the high electric conductiv- 
ity of the anode, the absolute value of electric potential is small. 
The electric potential in the porous cathode is a positive value 
and decreases from the cathode/electrolyte interface to the cath- 
ode current collector. Because the porous cathode is thin and its 
conductivity is relatively high, the variation of electric potential 
at cathode is very small. As shown in Fig. 6(c), the potential 
in electrolyte varies more significantly in the r-direction than 
in z-direction, and the current direction in the electrolyte is 
mainly vertical to the cathode/electrolyte interface. The reason 
is that the whole cathode is connected to the cathode current 
collector, and the potential difference between the electrolyte 
and cathode is much larger than that inside the electrolyte at 
z-direction. 

Fig. 7 illustrates the distribution of electromotive force 
(Nernst potential) along z-direction. ENernst decreases along the 
flow direction because of the Hz depletion and high water con- 
centration. The EMF is about 30 mV less at the fuel outlet than at 
inlet, implying that the reduction of the reversible cell potential 
is one of the major reasons for the cell performance reduction. 
Because the change of O2 molar fraction at porous cathode is 
relatively small, the primary effect to ENernst is H2 molar fraction 
and H20 molar fraction at anode/electrolyte interface. Compar- 
ing Fig. 4 with Fig. 7, the variation trend of Eyernst is consistent 
with that of H? molar fraction. 

Fig. 8 shows distribution of current density along the 
anode/electrolyte interface at z-direction. The current density 
along the cell z-direction is more uniform under the OC mode 
than under the other modes. The current density under the OC 
mode firstly decreases, reaches a minimum of 4200 A m~? in 
the middle part of a 2 cm cell and then increases toward the fuel 


(c) (d) 
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Fig. 6. Distributions of anode (a) and cathode (b) electric potentials, electrolyte (c) potential and the current direction (d) in the electrolyte. 
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Fig. 7. Distribution of ENemst along the anode/electrolyte interface for a 2cm 
long SOFC. 


outlet. The difference between the maximum and minimum cur- 
rent densities is only 400 Am~*, which is smallest among all 
these modes. The current density distribution is similar under 
the IC, BC and TC modes. For an example, the maximum of 
current density appears at inlet, and the minimum appears at 
outlet. Their current density at inlet under these modes are all 
close to 5700 Am~?. The minimum current density under the 
BC and TC modes is 4900 A m~ at the cell outlet, whereas the 
current density curve under the BC mode is concave because 
more energy is consumed by the longer current path. The cur- 
rent density under the IC mode at outlet is the smallest among 
all modes, and the difference of current density between inlet 
and outlet is 1500 Am~7. Because the key parameters affect- 
ing the electrochemical reaction rate are the local Hz and water 
concentrations, high H? concentration at the fuel inlet results 
in a high reaction rate. The current density under the OC mode 
keeps increasing when Hp concentration decreases in the fuel 
outlet region because the influence of the current path on the 
current density is more severely than Hz depleting rate. In a 
word, the electrochemical reaction under the OC mode is the 
most uniform one. 
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Fig. 8. Distribution of current density along the anode/electrolyte interface for 
a 2cm long SOFC. 
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Fig. 9. Distribution of anode overpotential along the anode/electrolyte interface 
for a 2 cm long SOFC. 


Fig. 9 and Fig. 10 show the distributions of anodic and 
cathodic overpotentials at electrode/electrolyte interfaces along 
the z-direction, respectively. The overpotential increases with 
current density increasing, vice versa. It is also seen that the 
cathodic overpotential is greater than the anodic overpotential 
in the cell because of inherently slower cathodic reaction kinetics 
[10]. 

Fig. 11 shows vector plots of current density inside anode 
and the effects of the four current collecting modes on the cur- 
rent path. The arrows express current-flow directions, and the 
arrow lengths are proportional to current density. Under the IC 
mode, current in anode flows from inlet to outlet and runs into 
cathode through the electrolyte, showing a gradually decreasing 
of current density in the outlet direction. Under the OC mode, 
current flow direction in anode is from outlet to inlet, and current 
density gradually decreases in current-flow direction. Under the 
BC mode, current in anode flows from both inlet and outlet to 
the tube center, and current density gradually decrease from the 
double ends to the tube middle. In this case, the current perpen- 
dicularly flows into the electrolyte in the tube middle. It can be 
seen that the longest current path under the BC mode is only 
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Fig. 10. Distribution of cathode overpotential along the anode/electrolyte inter- 
face for a 2 cm long SOFC. 
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(d) 


Fig. 11. Vector plots of current density inside anode under IC (a), OC (b), BC 
(c) and TC (d) modes. 


half of that under the IC and OC modes. For the TC mode, the 
current-flow direction is perpendicular to the anode/electrolyte 
interface, and the current path is shortest among all modes. As 
far as the current collection is concerned, the TC mode gives the 
maximum efficiency. The efficiency and loss of the different cur- 
rent collecting modes are calculated by the following efficiency 
loss equations: 


I 
efficiency = —. (14) 
TC 
I; 
loss = 1 — —. (15) 
Itc 


3.3. Effects of different current collecting modes 


Fig. 12 shows the cell efficiency loss under the IC, OC and BC 
modes and the ratio of the IC loss over BC loss as a function of 
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Fig. 12. The cell efficiency loss and IC loss/BC loss as a function of anode tube 
length. 


tube length for a 1.49 mm interior diameter micro-tubular SOFC 
at 800°C. It can be seen that with the tube length increasing, 
the efficiency loss increases due to the increasing anode current 
path. It is shown that the cell efficiency loss is the same under 
the IC and OC modes but much smaller under the BC mode. 
For the different length anode, the efficiency loss is about 2—6- 
fold lower under the BC mode than under the IC and OC modes. 
Since the efficiency loss increases sharply when the anode length 
increases, the value of IC loss/BC loss decreases. It can be con- 
cluded that the BC mode is the most effective actual current 
collecting mode, under which the efficiency loss is less than 5% 
on a 2.5 cm long cell. It is shown that a shorter cell under the BC 
mode can effectively reduce the efficiency loss. The selection of 
the anode tube length and current collecting mode is very crucial 
to minimize the efficiency loss. 


4. Conclusions 


In this paper, a 2D numerical analysis of the fluid dynamics 
and electrochemical kinetics has been conducted for anode- 
supported micro-tubular SOFCs. The variations of specie 
concentrations and electric potential in 2D have been obtained. 
The simulated current path inside the anode-supported micro- 
tubular SOFC could well reflect the current flow in an real 
cell. The current direction in electrolyte was mainly vertical 
to the cathode/electrolyte interface. The current flow in anode 
was mainly parallel to the anode/electrolyte interfaces and then 
flowed diagonally into the electrolyte. In addition, the cell 
geometry was also considered in the calculation of the cell per- 
formance. The current collecting efficiency was estimated with 
the simulation results of the cell under different current collec- 
tion modes. The efficiency loss under the BC mode was about 
2—6-fold lower than those obtained from the IC and OC modes, 
and the BC mode was the most effective current colleting mode. 
The anode current-collecting mode and tubular length were cru- 
cial for the design of micro-tubular SOFC. Those calculations 
using the model provided the optimum cell type and minimized 
efficiency loss. 


Acknowledgements 


The authors gratefully acknowledge financial supports from 
the Ministry of Science and Technology of China (No. 
2004CB719506 and 2006AA05Z147) and Natural Science 
Foundation of China (No. 20676132). 


References 


[1] S.C. Singhal, K. Kendall, High Temperature Solid Oxide Fuel Cells: Fun- 
damentals, Design and Applications, Elsevier, Kidlington, Oxford, 2003, 
p. 1. 

[2] A. Yan, M. Cheng, Y. Dong, W. Yang, V. Maragou, S. Song, P. Tsiakaras, 
Appl. Catal. B 66 (2006) 64-71. 

[3] H. Zhu, R.J. Kee, J. Power Sources 169 (2007) 315-326. 

[4] Y. Du, N.M. Sammes, J. Power Sources 136 (2004) 66-71. 

[5] J. Pusz, A. Smirnova, A. Mohammadi, N.M. Sammes, J. Power Sources 
163 (2007) 900-906. 

[6] N.M. Sammes, Y. Du, R. Bove, J. Power Sources 145 (2005) 428—434. 


254 D. Cui et al. / Journal of Power Sources 174 (2007) 246-254 


[7] H. Zhu, A.M. Colclasure, R.J. Kee, Y. Lin, S.A. Barnett, J. Power Sources 
161 (2006) 413-419. 

[8] T. Suzuki, T. Yamaguchi, Y. Fujishiro, M. Awano, J. Power Sources 163 
(2006) 737-742. 

[9] E. Achenbach, J. Power Sources 49 (1994) 333-348. 

10] M. Iwata, T. Hikosaka, M. Moriya, T. Iwanari, K. Ito, K. Onda, Y. Esaki, 
Y. Sakai, S. Nagata, Solid State Ionics 132 (2000) 297-308. 

11] H. Yakabe, T. Ogiwra, M. Hishinuma, I. Yasuda, J. Power Sources 102 
(2001) 144-154. 

12] H. Yakabe, T. Sakurai, Solid State Ionics 174 (2004) 295-302. 

13] K.D. Panopoulos, L.E. Fryda, J. Karl, S. Poulou, E. Kakaras, J. Power 
Sources 159 (2006) 586-594. 

14] R. Bove, S. Ubertini, J. Power Sources 159 (2005) 543-559. 


[15] S.H. Chan, K.A. Khor, Z.T. Xia, J. Power Sources 93 (2001) 130-140. 

[16] P. Costamagna, A. Selimovic, M.B. Del, G. Agnew, Chem. Eng. J. 102 
(2004) 61-69. 

[17] C.-Y. Chung, Y.-C. Chung, J. Power Sources 154 (2005) 35—41. 

[18] B.A. Haberman, J.B. Young, Int. J. Heat Mass Transf. 47 (2004) 
3617-3629. 

[19] Y. Ji, K. Yuan, J.N. Chung, Y.-C. Chen, J. Power Sources 161 (2006) 
380-391. 

[20] Y. Wang, F. Yoshiba, T. Watanabe, J. Power Sources 161 (2006) 380- 
391. 

[21] B. Todd, J.B. Young, J. Power Sources 110 (2002) 186-200. 

[22] R. Bove, P. Lunghia, N.M. Sammes, Int. J. Hydrogen Energy 30 (2005) 
181-187. 


